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^ ■ We present the results of a non-perturbative determination of the pion ma- 

trix element of the twist- 2 operator corresponding to the average momentum 
of non-singlet quark densities. The calculation is made within the Schrodinger 
functional scheme. We report the results of simulations done with the stan- 
dard Wilson action and with the non-perturbatively improved clover action 
and we show that their ratio correctly extrapolates, in the continuum limit, 
to a value compatible with the residual correction factor expected from per- 
turbation theory. 
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The calculation of hadronic matrix elements of Wilson operators entering the 
light cone expansion of two electroweak currents, directly connected to the moments 
of parton density distributions, requires non-perturbative tools. Lattice estimates 
have been produced so far for only the first few such moments. The results for the 
second moment of non-singlet quark densities, corresponding to the average momen- 
tum, are typically higher than the experimental values. However, these results have 
been obtained at fixed lattice spacing only and using, as renormalization factor for 
the bare operator, the one obtained from a perturbative calculation. It is the aim 
of our alternative approach to remove these approximations. 

The scope of the present paper is to provide values of the bare matrix ele- 
ment at different lattice spacing. After multiplication with the corresponding non- 
perturbative renormalization factor, these values can be used for an estimate of the 
continuum limit of the renormalized matrix element. In a series of papers, other 
essential ingredients needed for such an estimate were obtained. In particular, the 
Schrodinger functional (SF) renormalization scheme was discussed in [jlj, the non- 
perturbative step scaling function describing the scale evolution of the continuum 
renormalization factor in 0, the universality of the continuum limit of the step 
scaling function in || and the definition of a "renormalization group invariant" step 
scaling function in [f|] . The pion matrix element has already been calculated for the 
Wilson action and with standard periodic boundary conditions @: here we want 
to show how the SF can be used to reliably calculate not only the renormalization 
factor but also the physical matrix element itself. 

The use of the SF for extracting hadron correlation functions has been initiated 
with the calculation of pseudoscalar and vector masses and decay constants by the 
ALPHA collaboration ||. The calculation with the SF of a two-quark matrix ele- 
ment has never been attempted before and, with respect to a traditional method, 
presents two advantages. The first is the possibility of constructing "smeared" and 
gauge- invariant states as particle sources at the boundary (for the pion in this case). 
The second advantage is that the choice of SF boundary conditions allows inser- 
tion of the operator at T/2, where T is the total time extent of the lattice, to be 
compared with the maximum value T/4, reachable with ordinary periodic boundary 
conditions. If the operator of interest is inserted at a fixed physical distance, this 
amounts to needing a lattice for the SF that is a factor of 2 smaller in the time 
direction than the one using periodic boundary conditions. 

In our previous work we have evaluated a non-perturbative renormalization factor 
of the following operator: 

1 _ «- r 3 

12 (x) = -ip(x) 1{1 D 2} y^(ar) , (1) 

where _D M is the covariant derivative and the bracket around indices means sym- 
metrization. The calculation of the matrix element of this operator requires a non- 
zero pion momentum. In practice, this leads to a very noisy signal for large time 
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Figure 1: The correlation function 

separation. 

This has been observed also with periodic boundary conditions in ref. 0, but 
is even more crucial in our case, where we can reach a larger distance in time to 
project on the lightest state. We therefore decided to follow ref. || and compute 
numerically on the lattice the desired matrix element between pion states with a 
different lattice representation of the twist- 2 operator than the one defined in eq. (|]). 
This amounts to taking: 
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The advantage of this operator is that it can be computed at vanishing external 
momentum and is hence expected to show a much better signal to noise behaviour 
than 0i2 and not to be contaminated too strongly by lattice artefacts. Indeed, as in 
U , we find also within our SF calculation, the correlation function of O 00 to be much 
less noisy than 12 . In addition, we want to remark that we find with comparable 
statistics, the error of the matrixelement of the operator O o to be considerably 
smaller when using SF than for periodic boundary conditions as used in ||. 

The operator is inserted between two SF states defined at the two time bound- 
aries, as depicted in fig. |l[ We start by defining the boundary operators at the time 
boundaries and T: 



°o = ^EC(y)75yC(z) 



°r = h y E?w , y*T? 
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and then the correlation function (of mass dimension equal to V 



f M (x )=a i J2(°o0 o(x)0 T ). (4) 

X 

We also consider, for purposes of normalization, the correlation function 



h = -(0 o 0t) 



(5) 



Closely following the discussion in ref. ||, we easily see that for a time extent 
large enough, the matrix element is expected to be taken among the lightest particles 
coupled to the surface states, i.e. between pions, in our case. The general formulae, 
including the first excited state, are: 

fi „ p 2 e - m .T 

/m(x ) ^ p 2 e- m * T (n\O 00 \iT){l+p'e- Ax °+p"e- A ( T -^} . (6) 

The actual matrix element is obtained from /m(#o) only after a suitable normaliza- 
tion by /i, which takes out the effects (wave-function contribution) of the boundary 
quark fields. Assuming that there is a plateau region where /m(^o)//i — const = 
(ti\Oq \it) , and in which the first excited state gives essentially no contributions, we 
obtain the physical matrix element by a suitable normalization (see ||): 

(x) = — WOook) , (7) 

with k the hopping parameter appearing in the fermion action. The value of the pion 
mass is obtained, following ref. ||, from the time dependence of the pseudoscalar 
(/p) and/or axial-vector (/a) correlation functions. 

The set-up for our numerical simulation, performed in the quenched approxima- 
tion, was to choose lattices of physical size L 3 • T, with T taken to be T « 3 fm in 
all calculations. We performed two sets of simulations, one using the Wilson action 
and a second, using the non-perturbatively improved clover action. The correlation 
function projects on a single pion state at a distance in time corresponding to about 
1 fm. We then selected a time interval of 1 fm again to extract the matrix element. 
An example for the plateau behaviour of our correlation function and a fit to a 
plateau value with a distance of about 1 fm is shown in fig. 0. 

At each of our four values of fi = Q/g^, the calculation of the matrix element is 
performed at three values of the quark masses, using a multiple mass solver. Table |1] 
contains all our results for the matrix element (x) at various values of the quark mass 
and of j3, for both the Wilson and for the 0(a)-improved clover action. In fig. [^ we 
show one example for the chiral extrapolation of the matrix element. 

In table |^ we give the data for the matrix elements extrapolated to the chiral 
limit. They can now be used to test the continuum extrapolation, by building the 
ratio of the matrix elements calculated with the two different lattice actions we have 
used in our simulation. In fig. [| we show these ratios as a function of the lattice 
spacing. The data are extrapolated to the continuum as a linear function of the 
lattice spacing. Given the limited range in the bare coupling where the data were 
collected, the extrapolation procedure only extrapolates lattice artefacts that vanish 
like a power of the lattice spacing. Logarithmically dependent corrections, i. e. terms 
proportional to the bare coupling, cannot be taken into account in this way. Then, 
the result of the continuum extrapolation of the ratio is not expected to be one, 
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Figure 2: Example for the plateau behaviour of the correlation function 2K / /m 7r 
/m(zo)//i to extract the matrix element, taken at f3 — 6.2, k = 0.1346 on a 24 3 ■ 4£ 
lattice (clover action). 
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Figure 3: Example for the chiral extrapolation of the matrix element at (3 
using the clover action. 
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but the value of a correction factor evaluated at a bare coupling value of order 1, 
corresponding to the region in /3 where the extrapolation was made. The correction 
factor from a perturbative calculation at g 2 = 1 is 1.04 M. Figure [5] shows that 
our data are very well compatible with such a value and gives us confidence in the 
uniqueness of the continuum limit. 

We have also computed the renormalization factor Z 12 of the operator of eq. ([]]) 
at the smallest scale /i m i n we could reach in our previous work for the scale evolution, 
/Vin = l-436r with r ~ 0.5 fm. The values of Z\ 2 for our two lattice actions are 
shown in fig. | as a function of (3. In the figure the physical scale /x miri is kept 
fixed when varying (3 by suitably choosing the lattice size. In order to obtain the 
renormalization constant at exactly the values of (3 where the matrix element is 
computed, we performed an interpolation which is also given in the caption of fig. |j. 
We give the interpolated values of Z\ 2 in table 0. 

If we would now use Z\ 2 to renormalize the operator Oqq, the extrapolation to the 
continuum limit of the renormalized pion matrix element of the operator O o would 
acquire a similar correction factor, as in the case of the ratio of the bare matrix 
element as discussed above. This is why we do not proceed to extract directly from 
our matrix element its renormalized value in the continuum limit and we postpone 
to a forthcoming paper the non-perturbative evaluation of such a correction factor, 
which anyway is expected, from perturbation theory, to be of the order of a few per 
cent. 

We have demonstrated in this paper that the SF calculation of hadron matrix 
elements is feasible and could be applied to other interesting cases, such as the 
operators arising from effective weak hamiltonians. 
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Figure 4: Non-perturbative evaluation of Z 12 at a fixed inverse scale L = 1.436 r . 

The upper curve is for clover-improved fermions, the lower one for Wilson fermions. 

The numbers quoted in table ||] follow from the following interpolations (shown as 

dotted lines), which we take as a definition of the two renormalization constants in 

this range of (3 values: 

Zg = 1.2196 - 0.2244(/3 - 6) + 0.2117 (J3 - 6) 2 , 

Z$ = 1.0331 - 0.3570(/3 - 6) + 0.5714(/5 - 6) 2 . 



p 


K 


Action 


fit interval 


m^ 


(x) 


6.0 


0.153 
0.154 
0.155 


W 


11 -21 


0.4183(14) 
0.3594(15) 
0.2925(18) 


0.3096(17) 
0.3005(22) 
0.2903(32) 


6.1 


0.151605 
0.152500 
0.153313 


w 


12-24 


0.3663(14) 
0.3104(16) 

0.2528(19) 


0.3070(25) 
0.2967(33) 
0.2862(47) 


6.2 


0.150600 
0.151300 
0.151963 


w 


17-31 


0.3139(9) 
0.2669(10) 
0.2164(12) 


0.3035(19) 
0.2915(24) 
0.2769(34) 


6.3 


0.149259 
0.149978 
0.150604 


w 


24-40 


0.2977(10) 
0.2490(12) 
0.2007(15) 


0.3024(23) 
0.2874(33) 
0.2703(56) 


6.0 


0.1334 
0.1338 
0.1342 


c 


11 -21 


0.3988(11) 
0.3527(11) 
0.3005(13) 


0.2701(15) 
0.2658(18) 
0.2597(27) 


6.1 


0.1340 
0.1345 
0.1350 


c 


12-24 


0.3497(10) 
0.2919(12) 

0.2223(14) 


0.2729(16) 
0.2671(21) 
0.2612(40) 


6.2 


0.1346 
0.1349 
0.1352 


c 


17-31 


0.2796(7) 
0.2428(8) 
0.2009(9) 


0.2624(16) 
0.2569(20) 
0.2519(27) 


6.3 


0.1346 
0.1349 
0.1352 


c 


24-40 


0.2637(10) 
0.2281(11) 
0.1876(13) 


0.2677(26) 
0.2617(31) 

0.2552(40) 



Table 1: Our results for the two actions, Wilson (W) and non-perturbatively 
improved (C), that we have used. We also give the fit interval in xo/a corresponding 
to a distance of 1 fm. 
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Figure 5: The continuum extrapolation of the ratio of the bare lattice matrix 
elements computed with the clover-improved and Wilson fermion action. 



p 


a (fm) 


T/a 


L/a 


Z\2 


(x) 


N 


6.0 


0.093 


32 


16 


1.033 


0.2727(46) 


600 


6.1 


0.079 


36 


16 


1.004 


0.2683(67) 


573 


6.2 


0.068 


48 


24 


0.985 


0.2549(49) 


416 


6.3 


0.059 


64 


24 


0.977 


0.2474(76) 


391 


6.0 


0.093 


32 


16 


1.2196 


0.2471(44) 


600 


6.1 


0.079 


36 


16 


1.1993 


0.2535(48) 


800 


6.2 


0.068 


48 


24 


1.1832 


0.2405(41) 


500 


6.3 


0.059 


64 


24 


1.1713 


0.2427(61) 


319 



Table 2: Our results for the matrix element in the chiral limit at various values of 
the lattice spacing. We also give the renormalization factor Z 12 as obtained from 
interpolating the simulation data shown in fig. f|. The upper set of data belong 
to simulations with the Wilson action, the lower set to simulations with the clover 
action. In the last column we give the number N of configurations as used for 
evaluating (x). 
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